Measuring the Decoherence of a Quantronium Qubit with the Cavity Bifurcation 

Amplifier 
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Dispersive readouts for superconducting qubits have the advantage of speed and minimal invasive- 
ness. We have developed such an amplifier, the Cavity Bifurcation Amplifier (CBA) jjij, and applied 
it to the readout of the quantronium qubit Q. It consists of a Josephson junction embedded in a 
microwave on-chip resonator. In contrast with the Josephson bifurcation amplifier [13], which has 
an on-chip capacitor shunting a junction, the resonator is based on a simple coplanar waveguide 
imposing a pre-determined frequency and whose other RF characteristics like the quality factor are 
easily controlled and optimized. Under proper microwave irradiation conditions, the CBA has two 
metastable states. Which state is adopted by the CBA depends on the state of a quantronium qubit 
coupled to the CBA's junction. Due to the MHz repetition rate and large signal to noise ratio we 
can show directly that the coherence is limited by 1/f gate charge noise when biased at the "sweet 
spot" - a point insensitive to first order gate charge fluctuations. This architecture lends itself to 
scalable quantum computing using a multi-resonator chip with multiplexed readouts. 

PACS numbers: 



INTRODUCTION 

Quantum circuits based on Josephson junctions are 
candidates for the fundamental building block of a quan- 
tum computer - a quantum bit, or qubit [if. Several 
implementations have been tested @, 0, 0,M B S H]j 
which may be distinguished by the variable controlling 
the state of the qubit (charge, flux or phase) and the 
technique used for the readout. These systems can be in- 
dividually addressed, controlled and read, making them 
some of the most advanced solid state qubits. 

An enabling characteristic of many superconducting 
qubits is the existence of an optimal bias point where 
the qubit is immune to first order to fluctuations in ex- 
ternal control parameters. Nonetheless, even systems op- 
erated at this "sweet spot" have coherence times limited 
by second order fluctuations of the external control pa- 
rameters (llj . In this paper we describe measurements 
of a quantronium qubit operated at this optimal bias 
point 0, Q employing the newly developed Cavity Bi- 
furcation Amplifier (CBA) - a fast, dispersive, scalable 
readout system based on a driven non-linear supercon- 
ducting resonator [tj [l(| • With this architecture we can 
measure the fluctuations in the quantronium's coherence 
on time scales as short as a second, allowing us to probe 
the 1/f property of charge noise on these scales. These 
measurements confirm previous studies of the limitations 
of coherence times in charge qubits llj, [T2J and casts ad- 
ditional light on the fluctuating character of Ti itself and 
how its value depends on the measurement protocol. 

The non-linearity of the CBA is provided by a Joseph- 
son junction placed in the center of a A/2 on-chip copla- 
nar waveguide resonator with both an input and output 
coupling capacitor playing the role of Fabry-Perot cavity 




FIG. 1: (a) Scanning electron micrograph of a duplicate of 
the measured device, both of which are fabricated on the same 
chip simultaneously, (b) Schematic of our sample and its en- 
vironment. The sample consists of a split Cooper pair box 
whose loop is interrupted with a large junction (quantron- 
ium) . The state of the quantronium alters the switching prob- 
ability Poi of the CBA between its two metastable states, 
which can be inferred from the amplitude and phase change 
of transmitted microwave pulses through the device. 



mirrors (see Fig[TJ) When driven with a microwave signal 
at frequency Vd such that Vd < Uq (l — V3/Q), where Vq 
is the resonator small oscillation natural frequency and 
Q the loaded quality factor, this system can have two dy- 
namical metastable states which differ by their oscillation 
amplitude and phase. As the driving power P is ramped 
past the bifurcation power Pb , the CBA switches from 
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the state of low amplitude to the state with high ampli- 
tude. We detect the state of the non-linear oscillator by 
monitoring the amplitude and phase change of the mi- 
crowave signal transmitted by the resonator. In parallel 
with the CBA's junction we place a split Cooper pair box 
(SCPB), a circuit known as the quantronium, giving the 
resonator two bifurcation powers P^ and P^ < P^ 
depending on the state of the qubit, |0) or The two 
qubit states are mapped into the two metastable states 
of the CBA by ramping quickly the power P to a level 
intermediate between p] ^ and p] 1 ^ . If the quantronium 
qubit is in 1 1) , the CBA will switch to the high oscillating 
state, whereas if it is in |0) the CBA will remain in the 
low oscillating state. 

This readout has the advantage of being non- 
dissipative as the readout junction never switches into 
the normal state, unlike the original DC-biased quantro- 
nium readout [2J . This dispersive readout minimally dis- 
turbs the qubit state and since after switching we do not 
need to wait for quasiparticles to relax, the repetition 
rate is only limited by the relaxation time, Ti, of our 
qubit and the Q of our resonator. Like the DC readout, 
the CBA readout can latch [l7| , allowing enough time for 
the measurement of the complex amplitude of the trans- 
mitted wave, and therefore excellent signal to noise ratio. 
These characteristics were also present in the Josephson 
bifurcation amplifier 0, [uj , Il7| . which implemented a 
bifurcating non-linear oscillator using a lumped clement 
capacitor in parallel with the junction. However this ca- 
pacitor was fabricated using a CU/S13N4/AI multilayer 
structure which was difficult to fabricate and integrate 
with more then one qubit. Also the parallel plate ge- 
ometry suffered from inherent stray inductive elements. 
In contrast, the CBA is fabricated using a simple copla- 
nar waveguide geometry with no stray elements. The 
resonance frequency vq and the quality factor Q are con- 
trolled by the resonator length and output capacitor re- 
spectively. The CBA geometry thus offers the possibil- 
ity of designing a multi-resonator chip with multiplexed 
readouts, which would accommodate up to 10 qubits at 
once, an important step towards scalable quantum com- 
puting. The present work, in addition to the assessment 
of 1/f noise in a new architecture, is a first step in this 
direction. 



SAMPLE FABRICATION AND 
CHARACTERIZATION 

The resonator is initially fabricated using photolithog- 
raphy on a bare Si wafer 14|. A LOR5A/S1813 optical 
resist bilayer is used and the development is optimized 
to have at least 50nm of undercut beneath the S1813 to 
avoid wavy edges and to obtain a sloped edge on the 
resonator. This sloped edge is obtained by evaporat- 
ing a 200nm thick Al layer onto the sample at 0.2 nm/s 



with an angle of 5° and with a stage rotation of 10°/s. 
Next the quantronium is fabricated using electron beam 
lithography. We use a MMA/PMMA resist bilayer and 
the Dolan bridge double angle evaporation technique to 
fabricate our junctions [lij]. For this sample the split 
Cooper pair box is fabricated first inside the resonator, 
followed by the readout junction in a separate step using 
new e-beam resist. Using a hollow cathode Ar ion gun we 
obtain an ohmic contact between the two e-beam layers 
and the resonator. 

For the qubit sample described in this paper, we have 
been working with a linear regime resonance frequency 
i/ = 9.64 GHz and a Q of 160. An on-chip twin (FigCQi) 
of the SCPB of the quantronium had a normal state re- 
sistance of 15 kfi, with small junction areas of 0.05 ^m 2 . 
The usual gate line of the quantronium is absent in our 
design since we use the readout lines to access both the 
RF and DC gate lines, another advantage of this config- 
uration. 



QUBIT CHARACTERIZATION 

We performed gate charge and flux modulations while 
keeping the qubit in its ground state, to check that we 
have flux periodicity and 2e charge periodicity, as shown 
in Fig O The readout was first operated in the weakly 
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FIG. 2: Gate charge and flux modulations of our device. We 
operated in the weakly non-linear mode (P <C ) and mon- 
itored the phase (gray-scale) as we varied the applied gate 
charge iV 9 = C g V a /2e and flux $ ($ = h/2e). The large 
ellipsoidal contours can be interpreted as induced transitions 
between the energy levels of the qubit at multiples of the read- 
out frequency. The green fitted lines are transitions between 
the and 1 energy levels at the readout frequency of 9.64 GHz, 
while the orange fits are for transitions between the and 2 
energy levels at twice the readout frequency, 19.28 GHz 

non-linear mode (P <C p]^ 1 ) and we measured changes 
in the phase of the transmitted signal as the gate charge 
and flux were varied, keeping the frequency fixed at the 
maximal phase response point (see Fig [J). Apart from a 
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slow background modulation due to the changing suscep- 
tibility of the ground state, we observe sharp contrast on 
contours of ellipsoidal shape. These can be interpreted 
as contours of constant qubit transition frequency coin- 
ciding with the readout frequency or its double, an effect 
similar to that observed by Wallraff et al. [6[. Using 
the well known expressions for the energy levels of the 
quantronium [3], we can reproduce the shape of these 
contours, within the uncertainty due to the low frequency 
gate charge and flux noises, and extract a Jospehson en- 
ergy of the SCPB Ej of 15 GHz and charging energy 
E CP = (2e) 2 /2C E of 17 GHz, where C E is the sum of 
the capacitances of the junctions in the SCPB and the 
gate capacitance. 




FIG. 3: Spectroscopy peak as a function of gate charge, N B . 
Staying at zero flux, we measure Poi while sweeping v s and 
stepping N g . The theoretical fit of the resulting sinusoidal like 
dependance of the peak with N g is given by the red dashed 
line with the fit parameters Ej — 15.02 GHz and Ecp = 
17.00 GHz. The vertical lines with no N g dependance are the 
excitations between qubit energy levels enduced at multiples 
of the readout frequency, similar to those seen in Fig [5] 



To get a more precise measurement of Ej and Ecp, 
we performed spectroscopy on the qubit by applying a 
weakly exciting 1 fis long pulse, followed by a latching 
readout pulse. The switching probability Poi between 
the two metastable states of the CBA is measured as the 
spectroscopic frequency v s is swept for each gate charge 
step, at zero flux. Due to the crosstalk between the 
readout and the qubit, we have to ensure that we have 
zero leakage of spectroscopic power outside our pulse. 
This is achieved by gating the LO on the mixers shap- 
ing our pulses with a pulse shape a few ns longer then 
the spectroscopic pulse. As a function of frequency, we 
find a peak in switching probability, whose position varies 
with gate charge with the expected sinusoidal-like shape 
shown in Fig[3] The theoretical fit, shown in red, refines 
the previous determination of Ej and Eqp to the val- 
ues Ej = 15.02 GHz and E C p = 17.00 GHz. Zooming 
in to the double "sweet spot", N g = 0.5, $/$o = 0, 
where the qubit is immune to charge and flux noise 
to first order, we measure a Lorentzian spectroscopic 



peak of width Az/qi = 2.6 MHz and a Larmor fre- 
quency vqi = 14.35 GHz. This gives a dephasing time, 
T2 = l/irAfQi of 120 ns. However, large charge jumps 
move the biasing point off the "sweet spot" causing the 
linewidth to be widened. More accurate measurements 
of T2 will be obtained from Ramsey fringes where we can 
follow the variation of T2 with time. 



QUBIT MANIPULATION AND COHERENCE 
EXPERIMENTS 



Once the qubit parameters are known, we can perform 
experiments on the qubit to determine the qubit 's quality 
in terms of its energy relaxation time T\ and dephasing 
time T%. An essential part of these experiments is the 
need to control the state of the qubit precisely. This is 
achieved by applying a microwave pulse to the gate with 
rectangular envelope, of amplitude A, and time length r. 
In a frame rotating at the gate pulse carrier frequency, 
v s , the qubit state then precesses at a frequency, v v given 
by 0: 
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FIG. 4: (a) Rabi oscillations in the switching probability Poi, 
as a function of gate charge modulation AiV g and gate pulse 
time length r. AiV 9 is calculated from the Rabi pulse envelope 
voltage A reaching the sample through the attenuation in the 
input lines and is plotted in terms of Cooper pairs, AN g = 
C g A/2e. Oscillations in the switching probability Poi are 
seen with both AN g and r. (b) Fitted Rabi frequency VB.au 
vs ANg. As expected from a two- level system, vg.au scales 
linearly with AN g . 
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FIG. 5: (a) 3000 Ramsey oscillations as a function of free evolution time At. Each trace is 2.1 /is long with 3 ns per step. They 
each take 0.35 s to acquire. We can see visually the variation of T2 for the different Ramsey fringes by noticing the variation 
in contrast in the fringes near 2 fis. (b) Sample data fit. We average 5 of the acquired Ramsey traces shown in (a) and fit to a 
decaying sinusoid to extract the T2 which is then plotted in (c). For this case we have a coherence time of 842 ns and a Ramsey 
frequency of 26.9 MHz. (c) Distribution of T2 for 3000 of the Ramsey traces (600 fits). The black dashed line is the result of 
a simulation of the free evolution decay of the Ramsey fringes with 1/f noise fluctuations on the gate, S q (ui) = a 2 In the 
simulation we used 10 times more points compared to the data to obtain a smoother curve, (d) Corresponding distribution 
of Ramsey frequencies at four different flux biasing points. Each distribution has 3000 Ramsey traces. The blue histogram 
corresponds to the data in (a), (b) and (c). The Ramsey frequency is extracted from the position of the maximum of the power 
spectral density of each decaying sinusoid. The distributions are lopsided to higher frequencies as would be expected from 
fluctuations in the gate charge around our operating point at the "sweet spot". The dashed line is the expected distribution 
assuming the same 1/f charge noise as in (c). 



where the Rabi frequency [l8| VR a u = 
2EcpAN g (0 |2V| 1) /h involves the gate modulation 
ampitude in units of Cooper pairs AN g — C g A/2e and 
charge operator matrix element between the excited and 
ground states (0 \N\ 1). For A = we have free evolution 
at the Ramsey frequency VR amsev — \v m - v s \. 

We begin by measuring the Rabi oscillations. The 
pulse sequence protocol involves a gate pulse with fixed 
v s = vqi and varying amplitude, A, and time length r, 
followed by a latching readout pulse. The sequence is 
repeated 10 4 times to measure the switching probability. 
The oscillations of the switching probability as a function 
of t and A are plotted in FigUJ The extracted frequency 
VRaU scales linearly AN g (Fig 2b) as expected from a 
two-level system. From the position of the first maxi- 
mum of the Rabi oscillations we can calibrate the pulse 
time length needed for a 7r-pulse, which drives the qubit 
from the ground state to the excited state. Using this 
7r-pulse we measure the exponential decay of the popu- 
lation of the excited state (data not shown) and obtain 
the relaxation time T\ which, during the course of an ex- 



perimental run, varied between 1.4 /is and 1.8 /is, with a 
persistence time of a few seconds. These values are com- 
parable to the results of Vion et al. Q and Siddiqi et al. 

The maximum constrast measured is about 60%, lower 
then the maximum contrast of over 99.9%, calculated for 
the ideal case of a non-relaxing qubit given the measured 
parameters of the resonator. This loss can be attributed 
to three main sources. First the qubit relaxes before the 
readout takes place, because of its finite Ti, resulting in 
a 10% loss in contrast. Second the qubit relaxes to the 
ground state as the readout voltage approaches the bifur- 
cation voltage, resulting in a further 25% loss in contrast. 
Other measurements (not described here) have suggested 
that this loss in contrast could be from Stark shifting the 
qubit to lower frequencies during readout, where it can 
come in resonance with spurious transitions, possibly due 
to defects in the substrate or in the tunnel barrier. The 
remaining loss could be accounted for by the fact that 
the transition between the two oscillating states of the 
CBA is broadened by more than a factor of 5 from that 
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expected. 

To measure the coherence time, T2, we follow a differ- 
ent pulse protocol in which we apply two ir/2 pulses sep- 
arated by a free evolution period of length At, followed 
by a readout measurement. An example of the resulting 
Ramsey fringes is shown in Fig[5jD. By fitting to an expo- 
nentially decaying sinusoid we extract T%. The frequency 
of Ramsey fringes is well fitted by the absolute value of 
the detuning \v s — ^01 1, yielding a precise measurement 
of v m = 14.346 GHz. We follow the time evolution of T 2 
and Vqi by recording a Ramsey fringe every 0.35 s (Fig 
UK). We observe a variation of T2 from 150 ns to 850 ns. 
The distribution of T2 is an asymmetric bell shaped curve 
peaking around 600 ns, with a long tail extending down 
to 150 ns (Fig[5j:). If we would average over a 17.5 mn pe- 
riod (3000 of the above Ramsey traces) we would measure 
an average T2 which converges to 500 ns, similar to the 
first Saclay result Q obtained with a qubit with a similar 
Ej/Ecp- The T2 fluctuations are correlated with fluc- 
tuations in the Ramsey frequency, which only fluctuates 
towards higher frequencies, giving lopsided distributions, 
as shown on Fig [5jl. At the "sweet spot" where we are 
working, variations in gate charge necessarily increase the 
transition frequency whereas variations in flux decrease 
it. Variations in critical current would supposedly keep 
the distribution of frequencies more symmetric. We can 
therefore conclude that charge noise, not flux noise, is 
the dominant source of decoherence in our sample. Fur- 
thermore, if we suppose that the charge noise is Gaussian 
with a spectral density that has the usually invoked 1/f 
form [201 given by S q (uj) — a 2 /\ui\, we can check if our P 



data can be explained by this model. This was carried 
out by directly numerically simulating the correspond- 
ing variations in transition frequency and calculating the 
Ramsey signal in the conditions of the experiment. The 
distributions of both the extracted Ti and vp 
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shown by the dashed lines in Fig [5fc and d. We obtain 
good agreement with the data for a noise amplitude of 
a = 1.9.10 _3 e, agreeing with the range of previously 
measured values of this noise intensity parameter. To 
reduce sensitivity to this charge noise we can make the 
energy levels of the qubit almost insensitive to charge by 
increasing the areas of the junctions in the SCPB and 
hence increasing Ej/Eqp- An Ej/Ecp of 8 could give 
a T2 in the ms range and hence this device would be T± 
limited. 



CONCLUSION 

We have successfully implemented an improved version 
of the bifurcation amplifier based on an on-chip CPW 
resonator as a readout for the quantronium qubit. It of- 
fers ease of fabrication and a larger range of operating 
parameters (wn. Q) compared to the original JBA im- 
plementation [8j. Using this readout which captures in 



real time the fluctuations in qubit parameters, we have 
demonstrated that the main source of decoherence in our 
sample is charge noise. By using a larger Ej/Ecp, we 
could reduce the curvature with gate charge of the levels 
of the Cooper pair box and we should be able to re- 
duce the charge noise induced decoherence [21| . Further- 
more, the CBA geometry is particularly well adapted to 
the multiplexing of the simultaneous readout of several 
qubits. We have started in this direction by successfully 
measuring the bifurcation of 5 CBAs with only one input 
and one output line. This configuration offers a path for 
scaling up of superconducting circuits up to several tens 
of qubits. 
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